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Abstract 

We continue our investigation of the two-dimensional Abelian sandpile model in terms of a 
logarithmic conformal field theory with central charge c = — 2, by introducing two new boundary 
conditions. These have two unusual features: they carry an intrinsic orientation, and, more 
strangely, they cannot be imposed uniformly on a whole boundary (like the edge of a cylinder). 
They lead to seven new boundary condition changing fields, some of them being in highest 
weight representations (weights —1/8, and 3/8), some others belonging to indecomposable 
representations with rank 2 Jordan cells (lowest weights and 1). Their fusion algebra appears 
to be in full agreement with the fusion rules conjectured by Gaberdiel and Kausch. 

1 Introduction 

After early appearances [1, 2] and since its first formal study in [3], much attention has been 
directed to logarithmic conformal field theories (LCFT). This has led to significant progress, both 
for abstract models (see the reviews [4, 5]) and for lattice realizations [6, 7, 8, 9, 10, 11, 12, 13, 14, 
15, 16, 17, 18, 19, 20, 21, 22]. The known examples support the view that the critical lattice models 
described by logarithmic conformal theories have some intrinsic non-local features, in contrast to 
usual equilibrium models defined in terms of local Boltzmann weights with no global constraints. 
These non-local features on the lattice have an echo in the continuum in the form of peculiar 
properties of LCFTs, like the non-algebraic, logarithmic nature of their correlation functions, their 
non-unitarity, and non-diagonalizable Hamiltonians. 

Until very recently, only a few isolated lattice models were known to be describable by a LCFT. 
These include the dense polymer model [2, 23, 21], dilute polymers [7], percolation [24, 25, 26, 27, 7, 
17, 28], the Abelian sandpile model (ASM) [6, 8, 10, 11, 12, 13, 14, 15, 16, 19] and the close-packed 
dimer model [18]. The situation has considerably improved since an infinite family of integrable 
lattice models has been found [20, 22, 29], which correspond to logarithmic extensions of the usual 
minimal models. 

The simplest, and most studied LCFT has central charge c = —2. A local conformal field 
theory with this value of the central charge has been defined in [30, 31], and is known as the triplet 
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theory due to the existence of three conserved currents. This theory therefore possesses an extended 
algebra with respect to which it is rational. Despite the fact that the Virasoro representation theory 
for logarithmic theories has been worked out in some detail [32, 33, 34], no consistent, complete 
field theory with c = —2 other than the triplet theory has been constructed so far. 

On the lattice side, two models are known to have c = —2: the dense polymer model and the 
Abelian sandpile model (certain aspects of the dimer model can also be described by c = —2, see 
[35] ) . Although none of them is described by the triplet theory, the question remains open however 
as to whether they fall in the same universality class, and are described by the same conformal 
theory. 

An important question is therefore to investigate the precise structure of the logarithmic theories 
underlying these two lattice models. They are presumably more generic than the triplet theory 
which, because of its extended symmetry, is special. They are also likely to be quasi-rational but 
not rational. 

It is our purpose in this article to further investigate the nature of the theory describing the 
sandpile model. We do this by defining new boundary conditions, and by exploring their various 
properties with respect to each other (fusion rules) and with respect to other known observables. 

2 The sandpile model 

We start with a brief description of the model [36]; for further details and more results, the reader 
is refered to [37]. 

Every site i of a rectangular grid C is assigned a height variable hi, taking the four values 1, 2, 3 
and 4. A configuration C is the set of values {hi} for all sites. Given a configuration Ct at time 
t, the discrete dynamics makes it evolve to Ct+i, in the following way. First the height variable 
of a random site i is incremented by 1, hi — > hi + 1, making a new configuration C' t . Then if the 
(new) height hi in C[ is smaller or equal to 4, one simply sets Ct+i = C' t . If not, all sites j such that 
their height variables hj are greater than 4 repeatedly topple, which means that hj is decreased by 
4, whereas the height of all the nearest neighbours of j are increased by 1. This toppling process 
stops when all height variables are between 1 and 4; the so obtained configuration defines Ct+\. 

The addition of a single sand grain may trigger a potentially very large avalanche: the toppling 
of one site transfers sand to its neighbours, which in turn can topple and so on. When the dynamics 
is run for long periods, the sandpile builds up, and is eventually subjected to an avalanche spanning 
the entire system. This correlates the height variables over very large distances, and makes the 
system critical in the thermodynamic limit. 

The stochastic dynamics described above is Markovian and defines a time evolution in the space 
of probability distributions on the set of configurations. The asymptotic state of the sandpile is 
controlled by the unique invariant distribution P£ . It can be shown that P£ is uniform on the set 1Z 
of so-called recurrent configurations, and vanishes on the others. To be recurrent, a configuration 
must satisfy certain global conditions, and this is where non-locality enters. For what follows, it 
will be enough to know that the recurrent configurations are in one-to-one correspondence with the 
spanning trees on C. 
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More important is what happens on the boundaries. In the above description of the model, all 
boundary sites are dissipative: under toppling, since they have three or two nearest neighbours, 
one or two grains of sand fall off the pile. Such dissipative boundary sites are called open, and 
define the open boundary condition. The presence of dissipative sites is essential for the dynamics 
to be well-defined, since otherwise the toppling process might never stop. It is convenient to think 
of the sand grains falling off the pile as being transferred to a sink site, connected to all dissipative 
sites. 

Boundary sites can also be closed, i.e. non-dissipative. In this case, a toppling site loses as 
many sand grains as it has neighbours, generally three (two for corner sites). As a consequence, 
the height of a closed site takes only three (or two) values, usually chosen to be 1, 2 and 3 (1 and 
2). 

All boundary sites can be independently chosen to be open or closed, insofar as there is some- 
where at least one dissipative site. If this is so, the dynamics is well-defined, and so is the set TZ, 
which will depend on the boundary conditions chosen. The transfer of sand in a toppling can be 
encoded in the toppling matrix A, 

' Zi if i = j is either a bulk site or a closed boundary site, 

4 if i = j 'is an open boundary site, 
Aij = < (2.1) 
—1 if i,j are nearest neighbours, 

k otherwise, 

where Zi is the number of neighbours of i. The matrix A is simply the discrete Laplacian on C, 
subjected to specific boundary conditions. 

In the infinite volume limit, the invariant measure P£ becomes a conformal field theoretic mea- 
sure, and homogeneous open or closed boundaries are conformally invariant boundary conditions. 
In what follows, we will work on the upper-half plane (UHP), with an infinite-dimensional matrix 
A, denoted by A op resp. A cl if the boundary is fully open resp. closed. The difference A op — A cl 
is diagonal, with entries equal to 1 on boundary sites, and elsewhere. 

These two boundary conditions, open and closed, have been well studied and are well understood 
(see [8, 11, 19] and below). 



3 Spanning trees 

As recurrent configurations on C are in correspondence with spanning trees on £*, the grid C plus 
the sink site, Kirchhoff's theorem (see for instance [38]) implies that their number is given by 

|ft|=detA. (3.1) 

Spanning trees are acyclic configurations of arrows: at each site i of £, there is an outgoing 
arrow, pointing to any one of its An neighbours (if i is dissipative, the arrow can point to the sink 
site). A configuration of arrows defines a spanning tree if it contains no loop. By construction, the 
paths formed by the arrows all lead to the sink site, which is the root of the tree. If walked in the 
opposite direction, the paths start from the sink site (the root), enter the lattice through the open 
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boundary sites, which are the only ones connected to the sink, ramify and eventually invade the 
lattice. 

Even though the explicit mapping is very complicated and non-local, the spanning trees provide 
a description of recurrent configurations which is completely equivalent to that in terms of heights. 
The non-local conditions satisfied by the height values of recurrent configurations are encased in the 
property of arrow configurations to contain no loop, which is indeed a global constraint. However 
the spanning trees are often more convenient for actual calculations. In particular, they allow us 
to think of boundary conditions directly in terms of the arrows forming the spanning trees. 

The simplest possibility is to force boundary sites to have their arrows pointing in the same 
direction, either left or right. This defines two new boundary conditions, which we denote by — ► 
and <— . The more arrows we prescribe, and the stronger the wind blows in the trees ! The results 
below show that a wind blowing on the edge has long range effects deep in the forest and alters 
the general pattern of arrows far from the boundary. In what follows we examine the properties 
of the boundary fields which effect a change from an open or closed boundary condition to one of 
these two arrow boundary conditions, and in the limit where the finite grid goes to the upper-half 
plane (UHP). In order to review the way this is done, we start with the known case where an open 
boundary condition is changed to closed, and vice- versa [8]. 

Consider the sandpile model on the discrete UHP, with open boundary condition all along the 
boundary. It has toppling matrix A op . We modify the model slightly by closing one boundary 
site i. The new model has a slightly different toppling matrix which, from (2.1), is simply given 
by A ncw = A op + B, where the matrix B is infinite-dimensional, but has a single non-zero entry, 
namely Ba = — 1. The determinants of A op and A new give the number of recurrent configurations 
in both models, and are clearly infinite. However their ratio is finite, and is actually given by a 
1-by-l determinant (because B has rank 1), where G op = (A op ) _1 , 

det A new 2 

-±f— - = det (I + BG op ) = 1 - G°f = - ~ 0.637. (3.2) 

detA op 7r y ' 

The number of recurrent configurations has decreased by a numerical (unimportant) factor, but 
otherwise, the closing of site i has no noticeable effect at large distances. In the scaling limit, the 
two models are identical. 

Things change if we close a large portion / of an open boundary, say n consecutive sites. 
Then the new model has toppling matrix A new = A op + B, with B equal to minus the identity 
matrix on / and zero elsewhere. The ratio of determinants is now an n-by-n determinant. It decays 
exponentially with n, like e ~ 2nG / w where G is the Catalan constant, because the effective free energy 
of a closed site is smaller than that of an open site by an amount equal to 2G/7T. Subtracting this 
non-universal contribution from the free energy, a universal power law remains [8] 

2C det A new 2 r 

e ~" detA° p = ^ det(I + BG ° P ^ 1 ~ An ' n » L (3 - 3) 

For n large, the closing of / has long range effects. In the scaling limit, the ratio of determinants 
can be viewed as the correlator of two boundary changing fields (0 op,cl (O) cl ' op (n)), from which 
one concludes that the boundary condition changing field op ' cl = ^ c1 >°p has dimension — \. 
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The same procedure may be used for arrows. If we want to force the spanning trees to contain 
an arrow pointing from i to i + 1, two neighbour boundary sites, we modify two entries of the 
toppling matrix by setting An = 4 + 5 and A^j+i = —1 — 5. In the limit 5 — > oo, det A contains a 
term linear in 5, whose coefficient yields exactly the number of spanning trees having the prescribed 
arrow. Consequently, we may define A new = A op + B, where the two non-zero entries of B are 
Bu = 5 and B^j+i = —S. Similar to (3.2), the ratio 

1 det A ncw 1 3 4 

lim — —7- — = lim - det (I + BG op ) = - - - ~ 0.227 (3.4) 

s^oo 5 detA°P 5^oo 5 y 1 2 vr v ; 

gives the fraction of spanning trees which have an arrow between the sites i and i + Exactly the 
same formula holds if an arrow is prescribed on a closed boundary, the only difference being in the 
use of the inverse Laplacian subjected to a closed boundary condition on the real axis. One finds 

1 Hpt A ncw 1 1 

lim - j - . = lim - det (I + BG cl ) = - ~ 0.318. (3.5) 
s^oo 5 detA cl 5->oo 5 y ' vr y ' 

The same remarks as above can be made: it is only when a large number n of arrows are forced 
on the boundary that long range effects develop. Prescribing n consecutive arrows amounts to 
changing the entries of the toppling matrix through a perturbing matrix B, of dimension n + 1. 
The determinant ratio has a term proportional to 5 n , whose coefficient gives the fraction of trees 
having the n prescribed arrows. As in the case of a change from open to closed, the asymptotic 
value of these determinant ratios can be evaluated exactly. 



4 Effect of the wind 



Suppose that we prescribe n consecutive right arrows in an otherwise fully open or fully closed 
boundary (when there is only one string of arrows, their direction, left or right, makes no difference). 
The boundary itself is the set of points {(x, 1) € 1? : x G Z}; the first arrow points from (1, 1) 
to (2, 1), and the last n-th arrow from (n, 1) to (n + 1, 1). We label these n+1 sites by an integer 
1 < i < n+ 1. 

The defect matrix B corresponding to this situation is zero except on this set of sites, where it 
equals 

/5 -S ■ 
5 -S 

S -S 



B(5) = 



\ 



-6 
/ 



(4.1) 



\o o o 

The limit of 5 _n det[I + B(5)G] for S — > oo reduces to the principal (n + 1, n + 1) minor of 
B(1)G, and is explicitly given by 

1 dctA 11 - 



il™ ¥ det a = ¥ det[I + B{6)G] = det[GhJ ~ - r "■ 



(4.2) 



where G is G op or G cl depending on the boundary condition surrounding the arrows. They are 
both invariant under horizontal translations, so that the previous determinant has the Toeplitz 



5 



form det(<7j_j). Using the integral representation of the Green matrix on the plane and the method 
of images, one finds that the coefficients a m = Jq W ^e tmk a{k) are the Fourier coefficients of a 
function a, which explicitly reads 

a°P(k) = (l-e ik ) |2-cosA;-^(3-cosA:)(l-cosA:)}, (4.3) 

for the open boundary condition, and 

for the closed boundary condition. 

The asymptotic evaluation of Toeplitz determinants has been initiated in the celebrated Szego's 
theorem, and has since then been generalized in many directions. In the present case, the functions 
a op and cr cl both contain a so-called Fisher-Hartwig singularity [39], for which the asymptotic 
analysis requires a fairly recent result by Ehrhardt and Silbermann [40]. 

Let a(k) be a function on the unit circle with the following form, for Re a > — |, 

a(k) = (2-2 cos k) a T (k), (4.5) 

where r(k) is a smooth univalent function, nowhere vanishing nor divergent. Then the asymptotic 
value of the Toeplitz determinant formed with the Fourier coefficients of a is given by [40] 

det^O^^-EHn^-^e^ 10 ^) , n » 1. (4.6) 

In this expression, (logr)o is the zeroth Fourier coefficient of logr, and E[a] is a constant whose 
explicit value can be found in [40]. This remarkable theorem exactly fits the above situation since 
a op and cr ocl have the form (4.5) with (3 = |, on account of 

1 _ e ik = e iik - n)/2 V2-2COS/C. (4.7) 

The other result, discussed in the previous section, about the closing of n sites within an open 
boundary, follows from the same theorem, with (3 = 0, a simpler case proved earlier by Widom [41]. 
In case the boundary is open, from the function cr op , one has a = (3 = ^, and 

r op (k) = 2 - cos k - ^(3 - cos k)(l - cos k), (4.8) 

whose average value is equal to (logr op )o = — ^jr- We therefore obtain that, among all the recurrent 
configurations of the sandpile model on the open upper-half plane, the fraction of those which 
correspond to spanning trees having the n prescribed arrows on the boundary, is equal, in the large 
n limit, to 

(4.9) 

The exponential decay has the same origin as the one discussed in the previous section, namely the 
fact that a site with a prescribed arrow has a smaller free energy than when the direction of the 
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arrow is not constrained. Dividing out by this non-universal term, the ratio tends to a constant in 
the large n regime, 

e"-ra.„, n>> , (410) 

In the other case, when the boundary is closed, the function a cl has the form (4.5) with (3 = 
^, a = 0, and 

cl / 3 - cos k 1 1 - cos k 

The exponential decay, controlled by (logr cl )o = — ^p-, is weaker than in the open case, as expected. 
However the main difference is that a universal power law remains since the ratio behaves, for large 
n, like 

e ^n\^g^ =E[aC l ]n -l/^ n>>1 (412) 

Let us now interpret these two results in the context of a conformal theory. The two cases involve 
the change of boundary conditions at two points, separated by a distance n. In conformal theory, 
these changes are effected by the insertion of specific boundary fields. Therefore the exponents 
which control the algebraic decay of the ratios of determinants, and 1/4 respectively for the open 
and closed boundary condition, should equal the sum of the dimensions of the related boundary 
condition changing fields. We will make the assumption that the potential field dimensions are 
contained in the (infinite) Kac table of the (p,p') = (1,2) minimal model, which has c = —2, and 
given by 

^ = 2^1. (4 ,3) 

The values allowed by this formula are, in increasing order, h = — |, 0, |, 1, ^, 3, . . . (part of the 
Kac table is displayed in Section 6). As the relevant conformal field theory is logarithmic [19], these 
scaling dimensions can specify either primary fields or logarithmic partners of primary fields. 

For the open boundary condition, the sum of the dimensions of the two boundary fields is equal 
to 0, so that the dimension of both fields, <p° v ^ and 0^> op , must be 0. It is consistent to set, in 
the large n limit, 

^ n = (<rno)<r ,op (n)) = E[a o P] _ (4 M) 

One should note that, since the arrow boundary condition is oriented, the two fields <jf v ^ = (p^' op 
(outgoing arrow) and </> op,< ~ = <j)^' op (incoming arrow) are different. 

In the closed boundary case, the same argument leads to fields cl, ~^ and 0^ ,cl with dimensions 
adding up to \. From the values of h allowed by the Kac table, the only possibility is that one 
field has dimension — | and the other has dimension |. This may look at first sight problematic 
because it prevents the identification of the ratio (4.12) with a 2-point correlator (0 cl ' _> (O)^) _> ' cl (n)), 
as was done in the open case, since it vanishes identically. However one remembers that for the 
dynamics of the model to be well-defined, dissipation is essential. In the present case, because 
the whole boundary is non-dissipative (and there is no dissipation in the bulk either), one has to 
introduce dissipation at infinity by hand [11, 19]. It has been shown that the field corresponding to 
the insertion of an isolated dissipation, called u>n, is a dimension field, logarithmic partner of the 
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identity. Therefore the correct conformal interpretation of the result (4.12) is that it corresponds 
to a 3-point function, 

e^ n 'rc d ^ (n)l = <^(0)^ d (n) Wjv (oo)> = E[*«]n- 1 /\ (4.15) 

At this stage, by analyzing the change from open or closed to outgoing or incoming arrows, 
we have identified four new boundary condition changing fields: two fields 4> op ^ and M ,op with 
dimension 0, and two fields cl,_> and (f)~*' cl , one with dimension — |, the other with dimension 
|, although we do not know yet which one is which. In addition to these four fields, there must 
be other fields which make the change from one arrow boundary condition to the other one, by 
reversing the direction of the arrows. We first determine their conformal weights. 



5 Conformal weights 

In this section, we consider all changes of boundary conditions among open, closed, left arrows 
and right arrows, and determine the conformal weights of the corresponding fields. Under the 
assumption that these fields are primary, the conformal weights can be obtained rather easily by 
analyzing appropriate 3-point functions. It will be enough to consider an open boundary, in which 
we change the boundary condition on two consecutive intervals, say on [0, L] and [L, L + n]. 

The effect of inserting two segments of different boundary conditions can be computed as before. 
The change of boundary condition on the two intervals is implemented by an appropriate defect 
matrix B, from which one computes the number of recurrent configurations (or spanning trees) 
for such a boundary as det[A op + B]. Dividing that number by det A op yields the corresponding 
fraction of recurrent configurations as det [I + G op B], of dimension L + n + 1. Since three changes 
of boundary condition are involved, the way this determinant behaves in the scaling limit should 
be given by 3-point functions, up to the usual exponential factors related to the boundary free 
energies. As we are assuming that the fields are (quasi-)primary, the general form of a 3-point 
correlator is 

(MO) ML) ML + n)> = const. J^T^-H. ' ^ 

Because the determinants to be computed are no longer Toeplitz, the sort of theorem we used 
in the previous section for the insertion of a single segment cannot be applied, and we therefore 
resorted to numerical computations. The values of the determinants depend on L and n. In all 
cases discussed below, the numerical calculations have been carried out for four fixed values of L, 
namely L = 20, 30, 50 and 70, and for each value of L, n was varied between 1 and 150 (although 
the smallest values of n were usually discarded in the analysis). By fitting the numerical data to 
the general form (5.1), the values of the conformal weights can be extracted. 

We first consider a string of L arrows next to a segment of length n of closed sites. This yields 
two different cases, according to whether the arrows point to the left or to the right. In both cases, 
the exponential factor by which we multiply the raw data is equal to e (^+2n)G/7r^ gee (3 3) and 
(4.10). 
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For the arrows pointing to the right, we find that the numerical data are well reproduced by the 
above 3-correlator for the values of the scaling dimensions equal to hi = 0, h% = | and /13 = — g. 
As the relevant 3-point function is 



— g (0 op ^(O)^' cl (L)^°P(L + n)) = const. (5.2) 

we may conclude that the field cj) changing a closed boundary condition into ingoing arrows has 
dimension |. From this it follows that the field (p op '~* transforming the open boundary condition 
into outgoing arrows cannot be the identity. 

This first result, combined with those of the previous section, implies that the field chang- 
ing the closed condition into outgoing arrows must have dimension — |. This can be confirmed by 
looking at the previous situation with the arrows pointing to the left. In that case the numerical 
data are consistent with the general 3-correlator for h\ = 0, and /12 = /13 = — g. We then have 

g '--3- (^(0)^' cl (L)^°P(L + n)) = const. n l '\ (5.3) 

which confirms that cl '~* has indeed dimension — |. Again the field <ft op >^ of dimension cannot 
be the identity. 

The last boundary condition changing fields to be identified are those which reverse the direction 
of arrows. There are two different cases: a change from right to left arrows, or from left to right 
arrows. We start with the first case. 

Let us assume that the change from right arrows — ► to left arrows <— takes place at the point 
x. On the lattice, this corresponds to a string of right arrows, the last one pointing to x, and 
then a string of left arrows, the first one also pointing to x. The site x itself can either be open 
or closed. This makes a difference if one remembers that the arrows eventually lead to the root of 
the spanning tree through some dissipative sites. If the site x is non-dissipative (closed), then the 
chains of arrows pointing to x from either side, must go back into the interior of the lattice in order 
to find their way out to the root. But if x is dissipative (open), then the arrows can go out to the 
root directly through x. The global pattern of the arrows making up the tree is then completely 
different in the two cases, so that the presence of dissipation, even concentrated at a single site, 
has strong effects. Therefore we distinguish two fields (^° ,P ^ and <p^ c ^ according to whether the 
point at which the transition takes place is dissipative or not. 

The scaling dimensions of these two fields can be determined by computing the determinants 
corresponding to the situation where an open boundary contains a string of right arrows immedi- 
ately followed by a string of left arrows, of length L and n respectively and where the transition 
site is either open or closed. The exponential factor related to the free energies is now equal to 
e 4(L+n)G/7i^ rp^g va j ues Q f the determinants multiplied by these exponentials, for L, n in the same 
range as before, are consistent with the following correlators, 
L n 



op op op 



(0 op >^(O)0^'^(L)<^' op (L + n)} = const., (5.4) 



•+ ((/> op '^(0)^"^(L)^ op (L + n)) = const. k±H. (5.5) 



op d °P Ln 
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As <j) op '^ = (j)^' op has weight 0, we conclude that the fields rf>~~* and <fi^ c ^ have dimension 
and 1, respectively. 

For the second case, we consider the change of arrows from <— to — k The sites where the 
transition occurs are non-dissipative by construction (they have their arrow pointing along the 
boundary), so there is only one type of boundary condition changing field This time the 

numerical calculations are consistent with a 3-point function going to a constant for large L and 
n, namely 

-~~*~Z. ^- - p - (^(0)^(L)^°P(L + n)) = const. (5.6) 

It implies that the boundary condition changing field 4>^'~^ has dimension 0. 

Therefore, associated with the four boundary conditions open, closed, — ► and we have eight 
boundary condition changing fields: four of dimension 0, two of dimension — |, one of dimension | 
and another one of dimension 1. 

We know the conformal weights of the fields but not their precise nature. We have assumed that 
they are primary, and therefore generate a highest weight representation of the Virasoro algebra. 
But this is not a complete statement in a logarithmic conformal theory, as we need to answer the 
following questions: do these highest weight representations stand on their own, and if so, with 
which type of reducibility properties, or else are they part of bigger indecomposable representations, 
and if so, which ones ? 



6 Virasoro representations and fusion 

In the section, we make a short review of the results obtained originally by Gaberdiel and Kausch 
[32] and independently by Rohsiepe [33] regarding irreducible and certain indecomposable Virasoro 
representations for c = —2, and their fusion. These authors have considered a much larger class of 
models, namely the models (l,i) with central charges c = 13 — Q(t + i _1 ) for t > 1 integer, but in 
what follows, we focus exclusively on c = —2, (mostly) following the notations of [32]. 

The infinite Kac table is filled with conformal weights h ryS given by h r>3 = ( 2r ~ s ^ ~ 1 for r,s > 1. 
The lower left portion of the table is displayed below, with r running horizontally rightwards, and 
s running vertically upwards. Every scaling dimension in this table appears an infinite number of 
times. One may check that the first two rows s = 1,2, in yellow, contain every number exactly 
once, as does the first column with the exception of the weight 0, which appears twice. 

With every h in this list is associated a highest weight representation Mh, also called a Verma 
module, built on a highest weight state \h) satisfying Lq\H) = h\h), and annihilated by the positive 
modes of the Virasoro algebra, L n \h) = for all n > 0. It generates the representation by the free 
action of the negative modes, 

Mh = span{L_ ni L_ n2 • • • \h) : m,n 2 ,...>0}. (6.1) 

The algebraic structure of the Verma modules Mh has been completely determined in [42]. 
Each Verma module Mh contains, in addition to \h) itself, an infinite number of highest weight 
vectors, called singular vectors. They generate submodules, which form a linear chain of nested 
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subspaces. If h = h TiS , the module M r ,s = Mh rs contains a singular vector at level rs (value of 
Lq — h r>s ), which generates a submodule M r + S ,s- Correspondingly one defines the quasi-rational 
representation V r , s = -Mr,s/-M-r+s,s by setting to zero this singular vector at level rs. 

When r, s run over the Kac table, one obtains an infinite number of highest weight representa- 
tions V r ,s- They are not all distinct: V T)S = V r > s ' if an d only if h T)S = h r > s i and rs = r's', implying 
that the Kac table contains each V rs a finite number of times (which depends on r, s). For instance, 
the representations V r ,i appear only once, but the representations V r ,2 appear twice in the table 
since V Ti 2 = Vi 2r f° r every r > 1. 

The representations V r , s for all pairs (r, s) with the same value of h = h rtS have different 
reducibility properties. Indeed a value of h T)S is attained for an infinite number of pairs (r, s), but 
with increasing values of rs. Therefore, the representations V r s become more and more reducible 
as one goes deeper in the Kac table: if h TiS = h r ' iS i and r's' > rs, then V r ' tS ' is the quotient of 
M r ^ s i = M r>s by a singular vector lying at a higher level than for V TiS , and this implies the inclusion 
V r ,s C V r ' )S '. 

The only irreducible representations V r , s are those for s = 1,2 (and for r = 1, s even by the 
above equivalence). The corresponding Verma modules have nested submodules given by 

M r ,s >M r +s,s > M r +2s,s > M r+3StS ► • • • (r > 1, S = 1, 2), (6.2) 

implying that the singular vectors appear respectively at level rs, (2r + s)s, (3r + 3s) s, . . . The 
characters of these irreducible representations are 

xZ(q) = rt^rr- nS ^ r > 1, « = 1, 2. (6.3) 

The representations V r)S with r, s in the Kac table exhaust the set of highest weight representa- 
tions which are degenerate, that is, which contain a null vector at a finite level. With each of them 
is associated a primary field satisfying a degeneracy condition at level rs. This in turn implies that 
the correlators containing this primary field satisfy a partial differential equation of order rs. 

The previous discussion shows no essential difference with the usual treatment of minimal 
models. New features arise when one considers fusions. Indeed a general result of [32, 33] is that 
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the fusion of the highest weight representations V T)S closes on a bigger set of representations, which 
not only contains the V r)S themselves but also new, indecomposable representations 1Z T) \. We briefly 
recall the basic structure of the representations 1Z r ,\ following [32]. 

The indecomposable representations H r i can be seen as made up of two reducible highest weight 
representations Vi^r-i and Vi,2r+i ( r — !)• The two representations are graphically pictured by 
two vertical lines, tied together by the action of the Virasoro algebra (arrows). The dots represent 
certain specific states, and crosses mark singular null vectors. 

For r > 2, the two highest weight states have different dimension. 
The left representation Vi,2r-i is a reducible highest weight subrepre- 
sentation within lZ r> i, with highest weight vector £ r . The first singular 
vector (fr r lies at level r — 1, but is however not null in Vi2r— l- The 
second singular vector (p' r , at level 2r — 1, is null both in Vi^r-i an d in 
IZri, so that 4> r generates an irreducible subrepresentation V r ,i of !Z r ,i- 




The normalizations of 4> r , (p' r can be fixed by requiring 



r-l 



+ 



■)&, & = (L r _ 1 + 



= 0. 



(6.4) 



The right representation Vi,2r+i has highest weight state ijj r , with the 
same conformal weight as (f> r . The descendants of tp r , represented on the 
vertical line above tp r , are obtained in the usual way by the action of the 
(strictly) negative Virasoro modes on tp r , to which descendants of £ r are 
added in an appropriate way. Among them are p r and p' r , respectively the first and second singular 
vectors of V\ t 2 r +i, when this representation is taken alone. In TZ r ^i, positive Virasoro modes map 
p r to descendants of <fi r , and p' r to descendants of cj>' r . Thus only p' r is singular in 7?. rj i, and set to 
zero, p' r = 0. 

Whereas Vi^r-i is a subrepresentation, Vi^r+i is not, as it generates states of Vi^r-i when 
acted on by Lq and the positive Virasoro modes (in addition to states of V\ t 2r+i itself). This 
off-diagonal action, making the full representation TZ r \ indecomposable, can be defined by 



L[ = P r £ r , L p ljj r = , for p > 2. 



(6.5) 



They show that the full representation TZ ri \ can in fact be generated from the single (cyclic) state 
Vv> although itself not a highest weight vector. The parameter (3 r cannot be absorbed in the 
normalization of the various states, and so is an intrinsic parameter which labels inequivalent 
representations. 

The only states in TZ r> i which are highest weight states (which correspond to primary fields) are 
£ r and cj> r . They generate subrepresentations Vi^r-i and V r ,i, respectively reducible and irreducible. 
The quotient representation Tl r ,i/Vi,2r-i obtained by setting £ r to zero, is equal to Vi,2r+i- It ma y 
be noted that TZ Tj i also contains an indecomposable subrepresentation, generated by p r , whose 
structure is equal to the quotient of 7£ r +i,i by its 4> T +i and p r+ \ states. Quotient representations 
will play an important role in the subsequent sections. 
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For r = 1, the highest weight states of Vi 5 i and Vi 3 have both 
conformal dimension 0. The one on the left generates an irreducible 
subrepresentation: its first singular vector, at level 1, is set to zero, 
02 = i-i^i = 0, implying that is the identity field. The other 
representation on the right has also a lowest state lo\ with conformal 
weight 0. The Virasoro modes act on it in a way similar to the case 
r > 1, except that there is no state below 4>i to which loi can be mapped. 
The defining relations are thus 

Lqlui = 4>\ , LpUJ\ = , for p > 1. (6-6) 



Again the full representation TZi t \ is generated from the single state lo\. As for r > 1, p± (here at 
level 1) is not singular in IZn, but p[, at level 3, is singular in TZ r ^\ and null. The only highest 
weight state of is (pi, the corresponding quotient Tv^i/Vi,! being equal to Vi,3. The vector 
pi generates an indecomposable subrepresentation, which has the same structure as a quotient of 

7*2,1- 

According to [32], the set of representations V r , s and !Z T ,i for r > 1 and s = 1,2, is closed 
under fusion. Relying on explicit checks on the first levels, the authors of [32] were led to a general 
conjecture regarding the fusion rules. As far as we know, this conjecture has not been proved, but 
has been verified in a number of examples, most notably in [21], where reducible highest weight 
representations are also considered, and in [22]. The fusion rules read [32] 

V ri ,l*V r2 ,i = V m ,i, V ri ,i*V, 2>2 = V m , 2 , (6.7) 

m m 

V ri ,2*V r2>2 = U m>1 , V ri ,l*K r2>1 = Q) ft m ,i, (6.8) 
m m 

V ri , 2 *ftr 2 ,l = [V m _i, 2 © 2 V m , 2 © V m+ i, 2 ], (6.9) 
m 

*^r 2 ,l = ©2ft m ,i ©7^+i )2 ], (6.10) 

m 

where the summations are over m = \r\ — r 2 | + 1, |ri — r 2 | + 3, . . . , ri +r 2 — 1, and Vo, 2 = 7£o,i = 0- 
Two remarks should be made regarding these fusions 1 . First, each indecomposable representation 
1Z T) i arising in these fusion rules has a definite and constant value of (3 r . The first few values of (3 r 
have been computed in [32], for instance /3 2 = — 1, P3 = —18. Thus whenever the representation 
7^-2,1 appears above, it is meant to be the one with parameter /3 2 = —1- The second remark is 
that if a given representation arises in a fusion product, any of its quotient representations can also 
arise. 

Indecomposable representations of the above type have already appeared in the conformal 
description of the sandpile model, since it has been shown that the (boundary or bulk) dissipation 
field belongs to the lowest level of a 72-1,1 representation [11], and that the bulk height fields are 
found on the first level of a (non-chiral) 72 2i i representation [19]. 



1 l am grateful to Matthias Gaberdiel for stressing these. 
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7 Field identifications 



In this section, we derive more constraints on the nature of the boundary fields found in Section 
5. Our analysis and proposal will be based on two main assumptions: (i) the boundary condition 
changing fields are primary, and degenerate at an as low level as possible, and (ii) they belong to 
representations V T)S or 1Z T) i reviewed in the previous section, or to a quotient thereof. Assumption 
(i) is a consequence of the usual interpretation of the lowest lying state of the Hilbert space TC a ^ 
as corresponding to the field changing the boundary condition from a to (3. There is no deep 
reason for the second assumption, except that the V r , s arid !Z r ,i appear to form a natural supply 
of representations. Let us note that, as a consequence, primary fields with fractional dimension 
necessarily belong to irreducible highest weight representations V r ,2- In the next section, a few 
checks will be presented that support these assumptions. 

We will not dwell much on the field that swaps the open and closed boundary conditions. It 
has already been much discussed in [8, 11, 19], where considerable evidence shows that this field, 
called /i, is a primary field of scaling dimension — |, and corresponds to the highest weight state of 
an irreducible representation Vi 2, degenerate at level 2. 

Apart from the field /j, the two easy cases are the other two fields with dimension — | and |, 
which we will call respectively fx' and v. For them, our two assumptions above lead directly to the 
conclusion that they are the highest weight states of the irreducible representations Vi 5 2 and 1^2,2) 
degenerate respectively at level 2 and 4. 

All other fields have integral dimensions. To determine their properties, we mainly use the 
constraints coming from their composition law, expressed by the fusion algebra: for any boundary 
condition 7, the fusion 4> a, ' y * (fp^ must close on fields which interpolate between the boundary 
conditions a and /3, and should in particular contain the boundary condition changing field 4> a ^ . 

Then from the fusion 

we readily conclude that (f) op '^ must belong to a representation TL\ \. Being of dimension zero and 
different of the identity, it can only correspond to the uj\ field in 72.i i. However its primary partner 
(f>i, namely the identity, decouples completely from all correlation functions since the identity in 
fact does not interpolate between two different boundary conditions. Stated differently, conformal 
transformations of a general correlator (. . . (p a,op (zi) (j) op '^(z2) (j)^'^(zs) . . .) will generate an inho- 
mogeneous term (. . . (p a ' ov (zi) (ff^'^{zz) . . .) which vanishes identically because of the mismatch of 
boundary conditions on [21,23]. Thus the primary field <p\ is null, and (ff p '^ becomes the highest 
weight state of Vi 3 = 7£i 1/V1 1 (as we shall see in the next section, the singular state at level 1 in 
V1.3 is not null; this makes sure that op,_> is not the identity). We call this field a, leaving the 
notation uj for the logarithmic partner of the identity. 
Likewise the fusion 

implies that <j) op '^ should be identified with the £2 field in a 72-2,1 representation. It cannot descend 
to the quotient by cp2 = -^-1^2 since this would mean that £2 itself is the identity. We set </> op ' <_ = £2- 
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open 


closed 






open 


id. 




a = [0] € Vi, 3 


6 = [0] G K 2 ,i 


closed 


a = \— -1 G Vi o 


id. 


a' = \—-\ G Vi o 










id 


£ 2 = [0] G 7^2,1 (center op) 
£ 3 = [1] G 7^3,1 (center cl) 




(X = [0] € Vl, 3 


M'=[-|]eVi )2 


a' = [0] G Vi, 3 


id. 



Table 1: Summary of the fields (j) a ^ which make a change of boundary condition from a (row label) 
to (3 (column). The numbers in square brackets denote the scaling dimensions. 

Next case is 4>^'~^ which can be seen as arising in the fusion 

0«-,d * ^ci,- = ^ ★ ^ = Vi, 2 ★ Vi )2 = ^n- (7.3) 

The same arguments as for </> op '^ show that </>^'^ should be identified with the lo\ field of 72i i, 
projected down to the quotient Vi )3 . Let us note that this is consistent with the fact that <- ' - * 
should also appear in 4>^> op *(p op '^ , equal to V 1)3 *V 1)3 = Vi,i © Vi, 3 © Vi, 5 [21]. We set <p^>^ = a'. 

cl 

The field of dimension 1, must be in 

cl 

The only primary field with that dimension is the field £3 in 7£ 3j i. So we set <^~ >,< ~ = £3. 
The last field to consider is of dimension 0. The only fusion to contain it is 

0->,° P * ^o P ,^ = ^ ^ ^ = ^ ^ ^ = 2Ui i e 2 ^ 2 i 21Zs i ^ x 

In the representations appearing in this fusion, there are only two candidates for a primary field 
of weight 0: the field £2 hi a T^2,i, or the field u% in a 7^i,i quotiented to a Vi |3 . However another 
consistency condition is that the fusion (fi op, ~* * 0^° ,P ^ should close on the representation 7^2,1 
containing (/> op >^. This rules out the second possibility if one uses Vi |3 * Vi j3 = Vi t i © Vi j3 © Vi^ 
[21]. One concludes that (j>~^ '*~ = £ 2 must be the lowest lying state in a 7^2,1 representation. 

Table 1 summarizes the physical interpretation of these fields, along with their scaling dimension 
and the type of representation they belong to. We add two remarks. 

Remark 1 : The four dimension fields cr, c',^2 and ^2> although members of different repre- 
sentations, all satisfy the same degeneracy condition at level 3, since they are annihilated by 
(L 2 _ 1 -2L_ 2 )/--i. 

Remark 2 : From the way they appear in fusions, the two indecomposable representations 7^2,1 and 
7£ 3j i found above should have parameters 02 = —1 and /3 3 = —18. The two fields £ 2 = </> op '^ 
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and £2 = 0^°' P< ~ are primary with zero dimension and degenerate at level 3. Both have a singular 
descendant at level 1, which however cannot be null states since that would identify £ 2 an d £2 with 
the identity. So the degeneracy level of £ 2 j£2 i s as l° w as possible. In contrast, this is not the 

cl 

case of £3 = In a fully-hedged 72-3,1, the held £3 is degenerate at level 5, but possesses the 

singular descendant ^3 at level 2. If the latter is set to zero, then £3 would belong to a quotient 
representation IZ^.i/cfts with two peculiarities. First, the parameter (3% loses its meaning, namely 
it can be reabsorbed in the normalizations of the non-zero helds. And second, in the quotient, 
the state ^3 an d all its descendants no longer have partners, since all the states of which they 
were the logarithmic partners have been set to zero. So the quotient is an indecomposable non- 
logarithmic representation. Following [28], this type of representation can be called reducible but 
indecomposable representations of rank 1. 

One can now check that the physical interpretation we give to these helds is compatible with 
the way they should compose under fusion. For instance 

v = cl,< - = V 2i2 e <j) c1 ^ * </r^- = // * = Vi, 2 * 72 3 ,i = V 2 , 2 2V 3 , 2 © V 4)2 , (7.6) 

and, using the results of [21], 

6 = op '^ = n 2 ,i e <f) op ^ * ^°' p ^ = a * £ 2 = Vi )3 * n 2 ,i = 72.1,1 © n 2 ,i © 7e 3 ,i. (7.7) 

In other cases, the composition law is satisfied provided one takes a quotient (see the second 
remark at the end of the previous section). Examples have been already encountered above, and 
another simple instance is provided by 

^o P ,ci ^ 0ci,o P = ^ * ^ = y 12 * y 1)2 =tti i;L . (7.8) 

This fusion contains the identity in the double quotient of 72i,i by (f>\ and pi, since 72-i,i/{0i, Pi} = 
Vi,i. Despite the fact that 72i,i by itself contains the identity, the double quotient is necessary 
because the OPE /i(z)/i(w) closes on an irreducible representation, as shown in [11]. 
Another example is 

<P C1 ^ * <^ cl = n*v = Vi,2 * V 2 , 2 = 72 2 ,i, (7.9) 

also expected to contain the identity. However no quotient of 72-2,1 is equal to the irreducible 
Vi,i. The only possibility is to quotient 72-2,1 by the singular state 4>2 = L_i£,2, which forces the 
identification of £2 with the identity. This quotient 72-2, 1/^2 is of the type discussed in the second 
remark above, indecomposable but non-logarithmic. 

8 Higher correlators 

In this section, we analyze higher correlators, and bring further support to the field identifications 
made in the previous section. We have investigated about a dozen 4-point amplitudes on the upper- 
half plane, but we merely present here four representative examples. They would in principle 
provide a check on the parameters (32,03 for the representations 72. 2 ,i and 72-3,1, but for reasons 
explained below, no value could be extracted. 
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8.1 Amplitudes with 7?-2,i 

There are a few calculable amplitudes which allow to study the fusion of V12 * V2 2 = 72-2,1 • One of 
them corresponds to the following situation. We consider the sandpile model on the UHP with an 
open boundary, in which we close the sites on a first segment, then force a string of right arrows 
on another consecutive segment, and finally close a third consecutive segment. It involves four 
changes of boundary conditions, say at z\ < z 2 < z% < 24, and corresponds to the following 4-point 
function, 

---"---+ ^ 1 " 3 ^ ^""op--- ~ <M(*l)/A*2M*3)M*4))- I 8 - 1 ) 

From the degeneracy of \x at level 2, we obtain the general form of the above 4-point correlator, 



where x = 2122:34 . 

Z13 Z 24 



/ ( \ if \ i \ ( \\ ( Z 12Z14Z2± V/ 4 



X 

a + (3 log ■ 



1 



^13*23^34' 

When Z12 — > 0, the fusion of fi and // should close on the primary field a with no logarithmic 
singularity, implying (3 = 0. 

In order to compare with lattice results, we have taken the lengths of the two closed segments 
to be L = Z21 = Z43, and that of the arrow interval to be n = Z32. In this situation, the variable x 
is equal to (j^) 2 = t 2 , in terms of which the 4-point function reads 

t* (L) v(L + n) n(2L + n)) = — - — J = «(— ^=J . (8.3) 

In actual calculations, we have fixed L = 20, 30, 50 and 70, and for each value of L, n has 
been varied between 1 and 150. The corresponding ratios of determinants have been numerically 
computed, and then divided by the appropriate exponential factors, equal here to e - 4 ( i + n ) G '/ 7r i The 
previous conformal result suggests that the data, if plotted against the variable t = j^p^, should 
collapse on a single curve. The data collapse is manifest in Figure 1, where the quantity 

rW A new 

e 4(L+n)GA detA ( 8 4) 

detA°P v ' ' 

is plotted as colour dots, in terms of the variable t. The solid black curve represents the function 
of t in the r.h.s. of (8.3), and for a value of the coefficient a fitted to the data (a ~ 1.068). This 
strongly supports the conclusion that indeed 

(lJ,{zi) n'{z2) v(zz) fi(z&)) = a f Zl2Zl4Z2i \ 1 _ (8.5) 

v 213^23^34 y 

When Z3 — > 24, it can be reduced to a sum of 3-point functions upon using the OPE of \x and 
v. Conformal invariance fixes its form to be 

v(z)fi(0) = cJ 2 ' 1 z- 1 / 4 {6(0) + -^logz^ 2 (0) + -^V'2(0) + |^logz(L_ 1 ^)(0) 

l. lp2 *P2 op 2 

3/3 2 + 8 2/r , ., n . z 2 , r . . . 3z 2 
-z ' ' 



16/3 2 



(L-i<h)(0) + (L- 2 6)(0) + {L-i^m + ■■■} (8.6) 
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Figure 1: Comparison of numerical results against the conformal prediction for the situations 
depicted in (8.1). The colour dots represent the lattice numerical results for (8.4), for n from 1 up 
to 150 and L = 20 (red), 30 (orange), 50 (green) and 70 (blue). The abscissa is t = jk^. 

where fa is the parameter associated with the representation 7^2, i- 

Expanding the 4-point amplitude in powers of Z34, and using the previous OPE leads to the 
following 3-point functions 

C^ 2 ' 1 (ti(zi)fj,' (.22)6(34)) = az^ , C^f' 1 (/x(^i)/i' '(^2)^2(24)) = 0, (8.7) 

CS' 1 (mC^'C^C*)) = ^4 /4 (— + — )• (8-8) 

The first one has the correct form for primary fields of dimension — |, — | and 0. The second 
and third ones are compatible with the first one and the relations 4>2 = £—16 an d -^1^2 = ^26 ■ 
Unfortunately the vanishing of the (//// <fo) prevents us to check that /?2 has the expected value 
— 1 . We note that the expansion of the third relation around Z\ = Z2 implies that a and its level 1 
singular descendant L_icr both have a non-zero 2-point function with ip2- This shows that L-\a is 
not null, as announced in the previous section. 

If one takes the limit Z2 — > Z3, the expansion in powers of Z23 yields similar results. However in 
this case, the primary field £ is expected to be the identity, so that its descendant <f> = clearly 
decouples since <f> becomes null. 

Others situations allowing to study the OPE of Vi 2 and V2,2 can be examined in the same way, 
with similar results (no logarithmic term, decoupling of the <f>2 field). Among them, the following 
situation involves the boundary field related to having a height equal to 1 at an open boundary 
site. We consider again the UHP with an open boundary, in which we insert a segment of right 
arrows, close another segment of sites and ask that a given site in the remaining open boundary 
have a height variable equal to 1. In field theoretic terms, it corresponds to 

zi z 2 z 3 h = l 
-"op" d h £""qp~ ~ {^ 1 )v{z 2 )li{z 3 )h 1 (z A )). (8.9) 

The fields a, fi and v are the fields discussed earlier, while the field h\(z) is the boundary scaling 
field corresponding to the lattice observable 5{h z — 1) — P° p , i.e. the deviation of the probability of 
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Figure 2: Comparison of the numerical (dots) and the conformal (solid) results for the probability 
difference (8.14) for lengths z 2 i = z 32 = L, taken to be equal to 20 (red) and 30 (green), and a 
varying position Z4. The left figure shows the results when the height 1 is at a distance n to the 
left of the closed segment, so that Z43 = — n — 2L with n between 1 and 250. The right figure shows 
similar results when the height 1 is at a distance n to the right of the segment of arrows, so that 
Z43 = n, again running from 1 to 250. For convenience, both the numerical data and the conformal 
curves have been multiplied by a factor n 2 , so that the discrepancy between the dots and the solid 
lines is actually much smaller than what it appears to be on the vertical scale. 



having a height 1 at z from the value it would have if the boundary was fully open (and no other 
field inserted). The scaling fields hi(z) for all four height variables on an open or closed boundary 
are known exactly [12, 14]. In particular the scaling field for the height 1 on an open boundary is 
proportional to the boundary stress-energy tensor 2 [14] 

, / s /3 80 512\ m , , , , 

hi{z) = (--^ + 7 -AT{z). (8.10) 



7T 3TT 2 97T 3 / 

As T(z) = (L_ 2 )I(z), the 4-point amplitude (8.9) is computed simply by applying the appro- 
priate differential operator on the 3-point correlator (g(z\) v{z 2 ) ^(z 3 )) = Cz\^ z^^ 2 z 2 ^ 4 ■ The 
result is, with A\ the numerical factor in (8.10), 

A\z r i 
(cr(zi) u(z 2 ) n{z 3 ) h^)) = z 12 {z 2 4 + 3z u ) + z 23 z 2 4 ((r(zi) u(z 2 ) fi(z 3 )) . (8.11) 

° z 14 z 24 z 34 L J 

The expansion for z 2 — ► z 3 and the OPE (8.6) yields the relations 

(ffCsi)6teOM*4)) = °> (^(^1)^2(^)^1(24)) = 0, (8.12) 

^rC^ 1 (a( Zl )Mz3)hi(z 4 )} = A 1 C-^. (8.13) 

^p 2 ^^14 ^34 

Because a has zero weight and since h\ is a descendant of a zero weight field, the field ip 2 behaves 
like a primary field in the 3-point amplitude. 



2 Boundary height variables can be given simple expressions in terms of symplectic fermions [14]. However since 
then it has been proved in [19] that the symplectic fermions cannot account for all aspects of the sandpile model. 
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Note that the ratio 

n / m n pop (cr(zi) v(z 2 ) fi(z 3 ) h^zi)) Aiz 23 r i 

Pl(z 4 D.C.) - Pi = / / w n / u = ~o T^T z l2{Z 2 i + 3Z 34 ) + 223^24 (8.14) 

{a{z 1 )v{z 2 )^{z3,)) 8z 14 ^| 4 L J 

is the probability that the height at the open boundary site Z4 is equal to 1, given the boundary 
conditions pictured in (8.9), minus the same probability but with a fully open boundary (-P° p ~ 
0.104 [23]). This probability has been computed numerically for different positions of the various 
fields. The results are shown in Figure 2, and compare well with the conformal formula. 

8.2 Amplitude with 1Z 3A 

As a second example, in fact the only 4-point amplitude which involves the fusion to 72-3 1, we 
consider the sandpile model on the UHP with an open boundary, containing a first segment of right 
arrows, then a second interval of closed boundary sites, and finally a third segment of left arrows, 
as depicted here, 

------- l 1 I 2 cl 1 3 4 "5p--- ~ {a{z 1 )u{z 2 )u{z 3 )a{z 4 )). (8.15) 



From the null descendant of a at level 3, (L 2 _ 1 — 2L^ 2 )L^\a = 0, one can find the general form 
of the 4-point function as the solution of a third order differential equation, 



(a( Zl ) u(z 2 ) u(z 3 ) a(z 4 )} = z~^ [a Q= - ^) + (3 + 7 ( log + 2^) }• (8.16) 



The OPE f{z 2 )iy{z 3 ) should close on fields that interpolate between right arrows and left arrows 
(with a closed central site). If one postulates that, among these interpolating fields, the one with 
smallest dimension is the boundary condition changing field, then one can conclude that the OPE 
l/ (z 2 )i'(z3) cannot expand on fields with dimension lower than 1. This readily implies f3 = 7 = 0, 
and we are left with 

(<t(zi) u{z 2 ) v(z 3 ) a(z 4 )) = a z~ 3/4 — (8.17) 

V x 

To compare with the lattice results, we have set z 2 i = Z43 = L, taken as before equal to 20, 
30, 50 and 70, and z 32 = n, varying between 1 and 150. Then x is equal to (j^) 2 = t 2 and the 
4-point function reads 

((7(0) v(L) u{L + n) a(2L + „)) = J?L Lz£ m ( 8 .18) 

This function of t (i.e. omitting the factor n~ 3 / 4 ) has been plotted in Figure 3, along with the 
numerical values of the following lattice data 

3/4 2(4L+n)GA det (a in) 

6 detA°P • 18 yj 

Again the data collapse is clear and agrees very well with the result (8.18) for a fitted value of 
a ~ 1. 
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Figure 3: Comparison of numerical data against the conformal result for the situations pictured in 
(8.15). The colour dots represent the lattice numerical results for (8.19), for n from 1 up to 150 
and L = 20 (red) , 30 (orange) , 50 (green) and 70 (blue) . The horizontal coordinate is t = j^— . 

As we did in the previous subsection, the structure of 7^3,1 can be confirmed. Matching the 
general OPE (where (j) 3 = — 2L_ 2 )(, 3 is the first singular descendant of £3) 

u(zH0) = C^zV 4 {U0) + f (£-i6)(0) + f/{L^){Q) - Jt^iog^ 3 (0) 

-4^^3(0) + ...} (8.20) 

with the expansion of the 4-point function (8.17) yields 

^■>(si)6teM*4)> = a(— + —) , (8.21) 

v Zl3 Z34 / 

C^ 1 (a(z 1 )(L. 1 C 3 )(z 3 )a(z 4 )) = a (\ \-\ , (8.22) 

V z l 3 z 34 ' 

(a{ Zl )<f> 3 (z 3 )a(z 4 )) =0, (8.23) 

c^(a( Zl )M^H^)) = -^f ( - 4- + + -4r - 4) • (8- 24 ) 

10 \ Z 13 Z 13 Z 34 ^13 2 34 z 3 4 y 

These relations are compatible with the structure of 7^3,1 . Indeed the second relation follows from 
the first one, and the fourth relation is compatible with the second one and LiV> 3 = The 
third equation follows from the first one and cf) 3 = (L 2 ^ — 2L_ 2 )£ 3 , 

(a( Zl ) 4> 3 {z 3 ) a(z A )) = {L\ - 2C_ 2 )(a( Zl ) &(z 3 ) a(z A )) = a{L 2 _ x - 2£_ 2 )^ 4 - = 0. (8.25) 

^13^34 

Again the decoupling of 4> 3 from this amplitude does not allow us to extract the value of the 
parameter (3 3 . 

8.3 An example with dissipation 

As a last example, we examine a situation on the UHP which involves boundary condition changes 
on a non-dissipative boundary. In conformal terms, and because the boundary remains non- 
dissipative, one has to add a dissipation field ujn at infinity. This provides another opportunity to 
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check the consistency of the conformal picture, and allows further checks on the field £3 discussed 
just above. 

We consider the UHP with a closed boundary in which we insert a string of right arrows, 
immediately followed by a string of left arrows. Moreover the site at the junction between the two 
strings of arrows is taken to be closed. This situation is described in the scaling limit by, 

Zl Z 2 Z 3 OO 

— • — ^ ~ (// (21)6(22)^3)^(00)). (8.26) 

We first consider this 4-point correlator with the dissipation field located at a finite position 24, 
and then let Z4 — > 00. 

From the degeneracy of // at level two, the general 4-point function satisfies a second-order 
differential equation (even though lon is the logarithmic partner of the identity, this equation is ho- 
mogeneous because the 3-point function (//(zi) £3(22) ^'{23)) vanishes in the absence of dissipation), 
whose general solution reads 

5/4 1 

{^(z 1 )C 3 (z 2 )h\z 3 )lo n (z 4 )) = J^y/x(l-x)\a + l3(J-^ + axcsmy/r^)\. (8.27) 

Z21Z32 I- V V 1 - x J ) 

Since the configuration of arrows is symmetrical around the point Z2, one expects the corresponding 
amplitude to be invariant under the exchange z\ <-> Z3. As this transforms x into 1 — x, the above 
amplitude cannot be invariant unless (5 = 0, so that 

5/4 

£3(22) /A**) «»n(z4)) = a^— y/x(l-x). (8.28) 

Z21Z32 

In the limit za — > 00, in which x = Zl2Z3i t = —'\t becomes 

* ' Z13-224 213 ' 

(V'( Z1 ) £ 3 (z 2 ) M'(2 3 ) ^(00)) = -_=. (8.29) 

On the lattice, we have computed numerically the effect caused by the insertion of arrows (the 
ratio of determinants is finite). We have taken the lengths of the two strings of arrows to be L = z 2 i 
and n = Z32 respectively. As before the value of L is fixed to L = 20, 30, 50 and 70, and for each 
value of L, n ranges between 1 and 150. The corresponding ratio of determinants have been divided 
by the usual exponential factor, here equal to e _2 ( z '+ n ) G '/ 7r . i n addition, the previous conformal 
result suggests to first multiply the data by z^ = L 3//4 and then to plot them against the variable 
t = Figure 4 shows the comparison between the function of t in the conformal result above, 

and the numerical values of 

Ao+ A ncw 

L 3/4 e 2(L+n)GA , ( 8 . 3 0) 

det A cl v ; 

The plots show an excellent agreement, for a fitted value of the coefficient a ~ i. 

As a last remark, one should mention that (^3, the singular descendant of £3 at level 2, also 
decouples here. This readily follows from the fact that the function (8.28) is in the kernel of the 
second order operator implementing the level 2 degeneracy condition on the correlators, 

( f i'{z 1 )Mz2)^(z 3 )uj N (z 4 )) = (£ 2 „ x -2£_ 2 )( f i'(z 1 ) Uz2)v'(z 3 )uj n (z a )) = 0. (8.31) 
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Figure 4: Comparison of numerical results against the conformal result for the situation depicted 
in (8.26). The colour dots represent the lattice numerical results for (8.30), and for n from 1 up 
to 150 and L = 20 (red), 30 (orange), 50 (green) and 70 (blue). The solid curve represents the 
function of t in (8.29). The abscissa is t = 

9 Conclusions and perspective 

To summarize, we have defined, in terms of the spanning tree variables, two new boundary condi- 
tions in the two-dimensional Abelian sandpile model, and we have explored their nature in terms 
of a logarithmic CFT with central charge c = —2. Together with the well-known open and closed 
boundary conditions, they lead to seven new boundary condition changing fields (j) a ^ ' , with scaling 
dimensions in the set {— |, 0, |, 1}. We have examined many 3- and 4-point amplitudes within 
the conformal setting, and have found a full agreement with the corresponding lattice data. This 
brings further support to the consistency and relevance of the conformal description, and adds new 
entries in this description. 

It is worth stressing again the peculiarities and unusual features that these new boundary 
conditions have. 

1. The arrow boundary conditions carry an intrinsic orientation. To our knowledge, this is 
the first instance of oriented boundary conditions. A direct consequence of this is that the 
boundary condition changing fields which involve arrows have a vanishing two-point func- 
tion, or, formulated in another way, the metric on the space of states is off-diagonal in the 
representation basis. 

2. They cannot be uniformly imposed on a whole boundary. The reason for this is very clear 
in the sandpile model. The recurrent configurations are in one-to-one correspondence with 
spanning trees. Since a uniform arrow boundary condition on a boundary introduces a loop, 
it cannot be part of a spanning tree and therefore does not correspond to an allowed height 
configuration in the sandpile. In the conformal description however, this is a rather strange 
and new situation. The full implications of this from the general point of view of Boundary 
CFT need to be clarified. 
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3. All boundary condition changing fields are primary, but two of them, £ 2 = op '^ and £2 = 
(/> _>0 ' P<— , belong to an indecomposable representation 72-2,1 with rank 2 Jordan cells. The 
physical meaning of this and the physical interpretation of the lowest logarithmic partners 
(the ip 2 an d V4 fields) remain to be understood. 

cl 

4. A third boundary condition changing field, namely £3 = cp~* , ^~, is the lowest lying state of an 
indecomposable representation 72-3,1. In 72-3,1 this field is degenerate at level 5 but possesses 
a singular descendant at level 2. Even though this singular descendant has been seen to 
decouple in a number of amplitudes, we have no conclusive argument that it is actually null, 
and so we leave this question open. The same question as in point 3 regarding the physical 
interpretation of the field ^3 remains. 

Besides the questions raised above, there is clearly a number of other problems that need be 
answered before we can pretend to understand the sandpile model with boundaries. Among the 
most pressing and important ones, one can first mention the classification of all the observables 
which either preserve or interpolate between the four boundary conditions discussed here. Then one 
should also classify all conformally invariant boundary conditions present in the sandpile model, 
find their conformal description and determine the corresponding spectra of boundary fields. This 
is obviously much more challenging since an infinite number of boundary conditions are expected. 

Acknowledgments 

It is a great pleasure to thank Vyatcheslav Priezzhev for precious discussions and suggestions during 
the early stages of this work. I am also grateful to Matthias Gaberdiel and Jorgen Rasmussen for 
useful discussions, patient explanations and for valuable comments on the manuscript. This work is 
partially supported by the Belgian Internuniversity Attraction Poles Program P6/02. The author 
is a Research Associate of the Belgian National Fund for Scientific Research (FNRS). 

References 

[1] L. Rozansky and H. Saleur, Nucl. Phys. B 376, 461 (1992). 

[2] H. Saleur, Nucl. Phys. B 382, 486 (1992). 

[3] V. Gurarie, Nucl. Phys. B 410, 535 (1993). 

[4] M.A. Flohr, Int. J. Mod. Phys. A 18, 4497 (2003). 

[5] M.R. Gaberdiel, Int. J. Mod. Phys. A 18, 4593 (2003). 

[6] S. Mahieu and P. Ruelle, Phys. Rev. E 64, 066130 (2001). 

[7] V. Gurarie and A.W.W. Ludwig, J. Phys. A 35, L377 (2002); Conformal Field The- 
ory at central charge c=0 and Two- Dimensional Critical Systems with Quenched Disorder, 
hep-th/0409105. 



24 



[8] P. Ruelle, Phys. Lett. B 539, 172 (2002). 

[9] J. de Gier, B. Nienhuis, PA. Pearce and V. Rittenberg, J. Stat. Phys. 114, 1 (2004). 

[10] M. Jeng, Phys. Rev. E 69, 051302 (2004). 

[11] G. Piroux and P. Ruelle, J. Stat Mech. P10005 (2004). 

[12] M. Jeng, Phys. Rev. E 71, 036153 (2005). 

[13] M. Jeng, Phys. Rev. E 71, 016140 (2005). 

[14] G. Piroux and P. Ruelle, J. Phys. A: Math. Gen. 38, 1451 (2005). 

[15] G. Piroux and P. Ruelle, Phys. Lett. B607, 188 (2005). 

[16] S. Moghimi-Araghi, M.A. Rajabpour and S. Rouhani, Nucl. Phys. B 718, 362 (2005). 

[17] M.A. Flohr and A. Muller-Lohmann, J. Stat. Mech. P12006 (2005); J. Stat. Mech. P04002 
(2006). 

[18] N.Sh. Izmailian, V.B. Priezzhev, P. Ruelle and C.-K. Hu, Phys. Rev. Lett. 95, 260602 (2005). 

[19] M. Jeng, G. Piroux and P. Ruelle, J. Stat. Mech. P10015 (2006). 

[20] P.A. Pearce, J. Rasmussen and J.-B. Zuber, J. Stat. Mech. P11017 (2006). 

[21] P.A. Pearce and J. Rasmussen, J. Stat. Mech. P02015 (2007). 

[22] N. Read and H. Saleur, Nucl. Phys. B 777, 316 (2007). 

[23] E.V. Ivashkevich, J. Phys. A 32, 1691 (1999). 

[24] H. Saleur and B. Duplantier, Phys. Rev. Lett. 58, 22 (1987) 2325. 

[25] R. Langlands, P. Pouliot and Y. Saint- Aubin, Bull. Am. Math. Soc. 30, 1 (1994). 

[26] J.L. Cardy, J. Phys. A: Math. Gen. 25, L201 (1992). 

[27] G.M.T. Watts, J. Phys. A: Math. Gen. 29, L363 (1996). 

[28] J. Rasmussen and P.A. Pearce, Fusion Algebra of Critical Percolation, arXiv: 0706. 2716. 

[29] J. Rasmussen and P.A. Pearce, Fusion Algebras of Logarithmic Critical Models, 
arXiv: 0707. 3189. 

[30] M. R. Gaberdiel and H. G. Kausch, Phys. Lett. B 386, 131 (1996); Nucl. Phys. B 538, 631 
(1999). 

[31] M.R. Gaberdiel and I. Runkel, J. Phys. A: Math. Gen. 39, 14745 (2006). 

[32] M. R. Gaberdiel and H. G. Kausch, Nucl. Phys. B 477, 293 (1996). 

25 



[33] F. Rohsiepe, On reducible but indecomposable representations of the Virasoro algebra, 
hep-th/9611160. 

[34] H. Eberle and M. Flohr, J. Phys. A: Math. Gen. 39, 15245 (2006). 

[35] N.Sh. Izmailian, V.B. Priezzhev and P. Ruelle, SIGMA 3 (2007), 001. 

[36] P. Bak, C. Tang and K. Wiesenfeld, Phys. Rev. Lett. 59, 381 (1987). 

[37] D. Dhar, Physica A 369, 29 (2006). 

[38] V.B. Priezzhev, Sov. Phys. Usp. 28 (12), 1125 (1985). 

[39] M.E. Fisher and R.E. Hartwig, Adv. Chem. Phys. 15, 333 (1968). 

[40] T. Ehrhardt and B. Silbermann, J. Funct. Anal. 148, 229 (1997). 

[41] H. Widom, Amer. J. Math. 95, 333 (1973). 

[42] B.L. Feigin and D.B. Fuchs, Funct. Anal. Appl. 17, 241 (1983). 



26 



